On the resonant optical bistability condition. 
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Abstract 

We address a two level system in an environment interacting with the electromagnetic field in 
the dipole approximation. The resonant optical bistability induced by local-field effects is studied 
by considering the relationship between the population difference and the excitation field. The 
diversity of various systems is included by accounting for the system self-action via the surface part 
of the Green's dyadic in the general form. The bistability condition and the exact solution of the 
steady state optical Bloch equations at the absolute bistability threshold are derived analytically. 
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Because of its underlying nature and possible applications in the field of all-optical pro- 
cessing, the optical bistability (OB) was a subject of an intense experimental and theoret- 



optica 

J, a 



ical research. LU, [2| In the early studies the use of a saturable absorber to induce OB was 
suggested. |3l LZj The phenomenon was demonstrated experimentally for a cell of sodium va- 
por enclosed in a Fabry Perot interferometer and excited by a cw dyelaser.js] Later it was 
conjectured that the local-field corrections alone could give rise to the mirrorless OB. [9] 



10Ml4j and observed experimentally. 15[] 



This type of bistability was extensively studied 

The practical interest in all-optical devices faded to some extent as their solid state 
counterparts proved to perform better in terms of the switching speed and device density. 
However, the OB and optical hysteresis remain of considerable interest from the fundamental 



standpoint clS db clear manifestation of a nonlinear lie; 
of bistability mechanisms were discussed recently. 16 



it-matter interaction. Some new types 



19] Besides, the question of the OB and 



hysteresis has received a renewed attention in connection with novel hybrid OD nanoscopic 
systems, e. g., an artificial molecules comprising a semiconductor quantum dot (SQD) and 
metal nanoparticles (see Refs. I20 



-23 



and references therein). 
In this paper we address only the mirrorless OB induced by the local-field effects on a two- 
level system (TLS) interacting with the electromagnetic field in the dipole approximation. 
The local-filed correction leads to a self-action of the system, which results in a nonlinear 
relation between the applied field and the one acting upon the system. This type of the 
OB mechanism can be relevant for a large variety of systems: dense 3D assemblies of two- 



level atoms 



US 



aggregates, 



26 



optically dense thin films of TLS |24 l |25| or films of linear molecular 
hybrid metal-semiconductor systems, 2(J-l23| or a more general case of 
a TLS in an environment involving dielectric and conducting surfaces, such as, a stratified 
media, a microcavity or a nanostructure. 

OB can occur within a range of internal system parameters and external field intensities; 
identifying these ranges is therefore an important problems and its analytical solution is 
desirable. To the best of the author's knowledge, so far, it has been solved exactly only for 
two particular cases. Thus, Friedberg et al obtained analytically the bistability condition 
for the limiting case in which the active mechanism of the feedback is the nonlinear Lorentz 
shift of the resonance in a 3D gas of two-level atoms, which resulted from the near-field cor- 
rections. [14I Ignoring these corrections, Zakharov and Manykin derived the exact bistability 
criterion in the other limit in which the self-action is due to the radiated secondary field in 
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a thin 2D film of TLS. 25] The case when both these fields contribute into the self-action 



has been studied only numerically so far. 



29) 



We consider the simplest optical TLS coupled to an environment. The self-action of 
the TLS due to the environment can be described by the surface part of the Green's dyadic 
evaluated at the position of the TLS. Therefore, in this case the nonlinearity is characterized 
by a complex number which should satisfy some condition for the OB to occur. Below, we 
derive such condition analytically. 
Following Refs. 



14 



25 



29 



we address the OB condition by considering the relationship 
between the population difference and the excitation field within the framework of the Bloch 
equations for the 2x2 density matrix p mn (m,n = 0,1) of a TLS. In the rotating wave 



approximation these equations read: 

Z = -j(Z 



QP* + fi*Pi 



p = -(T + iA)P + QZ 



(la) 
(lb) 



where Z = pn — poo is the population difference between the excited and the ground state 
of the TLS, P is the amplitude of the off-diagonal density matrix element defined through 
p 10 = — (i/2)Pexp(— iut), 7 and T are the relaxation constants of the population and the 
dipole moment, respectively, A = ojq — uj is the detuning of the driving field frequency u from 
the TLS resonance uq, and Q = fiE/h is the total electric field E (in frequency units) acting 
upon the system, while /x being the TLS optical transition dipole. The dipole moment of 
the system is p = — i P/x, where P is its complex amplitude. 

To calculate the total field Q the Maxwell's equations are to be solved for a particular 
geometry. However, for a single dipole in an environment the total field Q can be represented 
in the following general form: 



tt = n + n 2 , Vt 2 = -iGP 



(2) 



where Qq is the renormalized incident field and fl 2 is the secondary dipole field. The latter 
field is due to the feedback of the environment and is related to the corresponding response 



function 



3CH32I] . It is determined by the surface part of the Green's dyadic evaluated at the 



position of the dipole. The near-zone component of the secondary dipole field, governed by 
Re(G), can originate from the near-field corrections, [l^] while the far- zone contribution, 



governed by Im(G), can be due to the radiation. 



In hybrid nano-systems comprising a 
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FIG. 1: The RHS of Eq. (J3ID for G = 6, A = 1 (solid line), G = 4, A = 1 (dashed line) and 
G = 8i, A = (dotted line) 

SQD, the feedback is provided by the secondary reflected field of the optical transition dipole 
moment. Thus, independently of its physical mechanism, the self-action is determined by 
the single complex valued parameter G. The total field Q acting upon the TLS depends 
therefore on its state. Together with the nonlinearity of the TLS that can give rise to a 
variety of effects, such as the OB. 

Under the steady-state conditions Eqs. (GQ) read: 

\U \ 2 z+i (r - dZf + (A + G R Z) 2 



jr z r 2 



(3a) 



p _ % ^0 /oux 

- (T-G l Z) + i(A + G R Z) ' 1 } 

where G R = Re(G) and G\ = Im(G). The above equation highlights a known result of cavity 
quantum electrodynamics: the real part of the Green's dyadic shift the frequency while the 



imaginary part renormalizes the decay rate. 30M32] Both effects can be active mechanisms 
of the OB. 

Eq. (I3aj) is of the third order in Z and can therefore have three real roots for some values 



of A, T, Gr, and G\ {V is used as a unit rate from now on). These three solutions are 
different when the right hand side (RHS) of Eq. fl3a|) has a minimum and a maximum (see 
the solid line in Fig. [[]), corresponding to two real roots of the RHS derivative, which satisfy 
the following equation 

2|Gf Z 3 + (|Gf + 2G R A-2G I )Z 2 -(l + A 2 ) = . (4) 

A threshold for bistability occurs when Eq. (jlj) has a double root (merged extrema - see the 
dashed or dotted lines in Fig. CO). In this case the solutions of Eq. (l3al) are: 

|G| 2 + 2(G R A-G I ) Z x 
Z1-Z2- 3|^p , Z,--- (5) 

Using the above root and the condition —1<Z<0 one obtains an important constraint 
for the detuning A: 

G R A-d< |Gf (6) 

The formal condition of the existence of the double root fl5]) is determined by the following 
equation on A: 

(|G| 2 + 2G R A-2Gi) 3 -27|G| 4 (l + A 2 ) = , (7) 

which is cubic in A and can also have three real solutions. It can be shown that only two 
of them, say Aj < A 2 , satisfy the constraint fl6]) and yield the upper and lower threshold 
values of the population difference Z and the external field Q . Therefore, Eq. (J4]) has three 
different real roots within a window of detunings A x < A < A 2 . The absolute bistability 
threshold occurs when the window shrinks into a point, i. e., when a degenerate root 
Ai = A2 of Eq. ([7]) appears. After some algebra, one can obtain the following alternative 
for the existence of such double root: 

G R + 2 {G\ - 4Gi - 8) Gl + Gf (d - 8) = , (8a) 

Gi + 1 = . (8b) 



Only two of the solutions of the biquadratic equation ( IHal) satisfy the constraint (|6]). 
Analyzing the roots of Eq. fll]) in the vicinity of the solution to Eqs. ( I8bl) and (JTj), one finds 
that the system can be bistable below the line given by G\ = — 1, while being critical (i. e., 
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FIG. 2: The stability phase map: outside the gray area the system can be bistable. Solid line 
represents the phase boundary given by Eq. ([9]). See text for more details. 

Zi t 2 — ~ 1) 011 the line for |Gr| < a/3. Finally, the absolute bistability threshold is: 

G R = 8 + 4G T - G : 2 + 8 a/1 + - l<d<8 (9a) 
-v 7 ^ < G R < a/3, Gi = -1 (9b) 

For completeness, we provide also the expression for the corresponding detuning: 

|G| 4 -4G|(G I + 3)-4G I 2 (G I -l) 

|G| 4 + 8G! [G? + G|(G I + 2)] 1 j 

Equation ([9]) provides the condition of existence of three real roots of Eq. f)3al) . We 
studied the stability of these by analyzing Lyapunov exponents of Eqs. ([T]) in the vicinity 
of a stationary solution and calculated the stability phase map numerically. Figure |2] shows 
such phase diagram in the space of the feedback parameter G. Three different real solutions 
of Eq. f )3a|) exist within the white area, two of them being stable. The line dividing the two 

6 



phases represents the analytical condition Eq. (Q which gives also the bistability threshold. 



The general criterion Eq. ([9 



25] 



recovers both reported exact results: G\ > 8 for Gr = 
and Gr > 4 for G\ = |14| (the corresponding threshold values of G are marked by full 
circles in Fig. [2]). Note the mirror symmetry of the bistability map with respect to the line 
Gr = 0, which reflects the invariance of Eq. (15a|) under the simultaneous change of signs of 
G R and A [29]. 

Eqs. (jSJ) and ( fit)]) give the full solution of the steady state problem at a bistability thresh- 
old, while Eq. fl9]) is the bistability condition, which constitute the main result of the paper. 

Finally, we note that novel hybrid plasmonic nano systems, such as a SQD and a metal 
nanoparticle complex (see Ref. |23| and references therein) or a SQD embedded in a stratified 
medium are excellent model systems with adjustable nonlinearity. The self-action field in 
these complexes is the secondary reflected field of the SQD optical transition dipole moment 
acting back upon the SQD. In the most general case such environment feedback results 
in a complex self-action parameter G which can be engineered by an appropriate choice 
of materials, geometry and/or external control parameters. Thus, if the SQD transition 
frequency is far from the plasmon resonance of the system then, typically, |Gr| 3> |Gi| 
and the dominant mechanism of the bistability is the nonlinear Lorentz shift 14] . If the 
excitation frequency is near the surface plasmon resonance then the feedback can be almost 
purely imaginary: |Gi| ^> |Gr| with Gt > 0. Algebraically, this case is equivalent to the one 
of the radiation induced self-action [25| . A very interesting case is G\ < which can hardly 
be realized in the traditional extended 3D or 2D systems of two-level atoms. Nanoscopic 
hybrids are much more promising from this point of view due to their greater diversity. In 
the latter case various instabilities such as auto-oscillations can be expected under steady- 
state excitation. However, the detailed study of these instabilities goes beyond the scope of 
the present contribution and will be analyzed elsewhere. 

Summarizing, we addressed the mirrorless resonant optical bistability of a two-level sys- 
tem in an environment by considering the relationship between the population difference 
and the intensity of the excitation field. The feedback of the environment was introduced in 
a general from via the surface part of the Green's dyadic. We derived the analytical bistabil- 
ity condition and obtained the exact steady state solution of the Maxwell-Bloch equations 
at the absolute bistability threshold. Our findings open the possibility to easily analyze 
diverse physical systems, determine experimental conditions (such as, the appropriate de- 
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tuning from the bare resonance and the intensity of the external field) necessary to observe 
the bistability and other nonlinear optical effects, as well as design and engineer new systems 
with desirable nonlinear optical properties. 
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